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Abstract. We study the effect of interfacial phenomena in two-dimensional perfect and random (or dis¬ 
ordered) g-state Potts models with continuous phase transitions, using, mainly, Monte Garlo techniques. 

In particular, for the total interfacial adsorption, the critical behavior, including corrections to scaling, 
are analyzed. The role of randomness is scrutinized. Results are discussed applying scaling arguments and 
invoking findings for bulk critical properties. In all studied cases, i.e., g = 3, 4, and g = 8, the spread of 
the interfacial adsorption profiles is observed to increase linearly with the lattice size at the bulk transition 
point. 

PACS. 75.10.Hk Glassical spin models - 05.50-l-q Lattice theory and statistics (Ising, Potts, etc.) - 05.10.Ln 
Monte Carlo method, statistical theory 


1 Introduction 

At the interface between coexisting phases various intrigu¬ 
ing phenomena may occur, like roughening and wetting m 
An interesting complication arises when one considers the 
possibility of more than two phases. A third phase may be 
formed at the interface between the two other phases. An 
experimental realization is the two-component fluid sys¬ 
tem in equilibrium with its vapor phase [5] . The situation 
may be mimicked in a simplified fashion via multi-state 
models by fixing distinct boundary states at the opposite 
sides of the system. The formation of the third phase, with 
an excess of the non-boundary states, in such models has 
been called “interfacial adsorption” Em- 

Various aspects of interfacial adsorption, including its 
novel critical behavior at the bulk phase transition and the 
interplay with wetting, have been investigated quite exten¬ 
sively. In particular, specific models, like Potts and Blume- 
Capel models, have been studied using numerical tech¬ 
niques, especially, Monte Carlo (MC) methods and density 
renormalization-group calculations |3I5I6I7I8I9I10I11I12] . 
Furthermore, scaling and analytic arguments have been 
presented |3l6l8l13l14lTm6] . In particular, critical expo¬ 
nents and scaling properties of the temperature and lattice 
size dependencies have been determined. The fundamental 
role of the type of the bulk transition has been clarified, 
with isotropic scaling holding at continuous and tricritical 
bulk transitions, while interfacial adsorption is described 
by anisotropic scaling at a bulk transition of first-order 
type. 

Despite the resulting interesting insights, open ques¬ 
tions still exist. In this article, we shall mainly consider 


two aspects which have been largely neglected before, the 
adsorption profiles and the effect of randomness on the in¬ 
terfacial adsorption in two-dimensional (2D) Potts mod¬ 
els m Concretely, we shall consider the g-state Potts 
model on the square lattice, where the ferromagnetic near¬ 
est neighbor couplings may take the values Ji or J 2 , oc¬ 
curring with the same probability and being distributed 
randomly. Of course, when the ratio r = J 2 /J 1 takes the 
value I, one deals with the pure (or perfect) case. In the 
dilute (or random) case, the position of the interface, as 
well as the extent of the intervening third phase of non¬ 
boundary states, may be strongly affected by the spatial 
distribution of the couplings. Because of the inhomogene¬ 
ity of a given bond realization, it seems, obviously, worth¬ 
while to monitor not only the total excess adsorption of 
the non-boundary states but also the local structure of the 
interfacial adsorption, especially the adsorption profiles. 

The bulk transition temperatures of these Potts mod¬ 
els are known exactly from self-duality for arbitrary values 
of the internal states q and disorder-strength ratios r [15] . 
Accordingly, analyses on the critical behavior of the inter¬ 
facial adsorption, based on extensive MC simulation data, 
as it is also done in the present paper, are significantly 
simplified. 

Bulk criticality of such random Potts models on the 
square lattice has attracted much interest, partly, because 
the transition is of continuous type for all values of q, 
while being, in the perfect case, of first order for q > 
Amm . Then, the analysis of the interfacial adsorption in 
these models may be simplified by the fact, that isotropic 
finite-size scaling is expected to hold at continuous tran¬ 
sitions mm- Static and dynamic bulk critical properties 
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of the random Potts models have been estimated, suppos¬ 
edly, rather accurately, using a variety of, predominantly, 
numerical methods m- Some of the these results will turn 
out to be very useful in our study on the interfacial ad¬ 
sorption. 

Attention should be drawn to related previous work on 
interfacial phenomena in dilute ferromagnetic Potts mod¬ 
els, in particular, considering hierarchical lattices, i.e. ap¬ 
plying the Migdal-Kadanoff real space renormalization to 
the square lattice m, or performing a preliminary MC 
study for the square lattice model [21] . 

The outline of the article is as follows: In the next Sec¬ 
tion [51 the model and the methods, especially, MC simu¬ 
lations of Metropolis and Wolff type, will be introduced, 
followed by the discussion of our main results in Section [31 
The summary. Section |T 1 will conclude the article. 


2 Model and methods 

We shall study the nearest-neighbor random-bond q-state 
Potts model on the square lattice described by the Hamil¬ 
tonian 

'H = — ( 1 ) 

{(id)) 

where the Potts variable at site i, Si, takes the values 
i,2,...,g[n]. The ferromagnetic random couplings Jij > 
0 between nearest neighbor sites i and j are either, with 
probability p, Ji or, with probability 1 — p, J2- In the case 
Ji > J2, one has either strong or weak bonds. Of course, 
r = J2/J1 = 1 denotes the pure (or perfect) 9-state Potts 
model. In our study, we shall consider Potts models with 
9 = 3 , 4 , and 9 = 8. 

In this article, we shall consider the system at its self¬ 
dual point, where both couplings occur with the same 
probability, p = 1 / 2 . Then the phase transition temper¬ 
ature, /cbTc/Ji, between the ordered ferromagnetic phase 
and the high-temperature disordered phase is known to 
follow from [ 15 ] 

^g(Jl/fcBTc) _ 0( g(rJi/fcBTc) _l'^ ( 2 ) 

In the random case, the phase transition is of continu¬ 
ous type for all values of 9, while, in the perfect case, the 
transition is continuous only for 9 < 4 , being of first order 
at larger number of Potts states 9. Exact values of the 
critical exponents are, so far, only known in the perfect 
case m- Numerical analyzes, in the dilute case suggest 
that the bulk critical exponents depend rather mildly on 

9 m- 

The degeneracy between the 9 equivalent Potts states 
may be lifted by appropriate boundary conditions. Here, 
to study interfacial adsorption, we shall employ special 
boundary conditions, distinguishing the cases [1 : 1] and 
[1 : 2 ]. For the case [ 1 : 1 ], the Potts variable is set, at all 
boundary sites, equal to 9 = 1, while for the case [1:2], 
the variable is set equal to 1 at one half of the boundary 
sites, and to 2 at the opposite half of the boundary sites. 



Fig. 1. (color online) Typical equilibrium MC configuration 
of the L = 100 9 = 8 random-bond (r = 1/10) Potts model 
at temperature T — 0.98Tc showing an interface for a par¬ 
ticular bond-disorder realization. Red color depicts the q = 1 
states, blue color the q = 2 states, whereas the non-boundary 
states (q > 3) adsorbed at the interface are shown blackened. 
Note that the fixed boundaries [1 : 2] are also included in this 
illustration. 

Then, the boundary condition [1 : 2 ] introduces an inter¬ 
face between the 1-rich domain (or phase) and the 2-rich 
domain (or phase). By examining typical MC equilibrium 
configurations, as shown in Figure |TJ it is seen that at the 
interface between the 1- and and 2-rich domains an excess 
of the non-boundary states is generated compared to the 
case in the absence of an interface. 

Then, the interfacial adsorption, W, measuring the 
surplus of non-boundary states induced by the interface 
between the 1-rich and 2-rich regions, is defined, for lat¬ 
tices with non-boundary sites, with L being the linear 
dimension of the lattice, by [ 3 ] 

^ ~ Yj [(^Si,n)[l:2] ~ (^Si,n)[l:l]] j ( 3 ) 

n i 

summing over all non-boundary sites i and over all non¬ 
boundary states n = 3 , 4 ,... ,9. The thermal average is 
taken. Obviously, W may be interpreted as the effective 
width of the domain of non-boundary states between the 
1- and 2-rich domains. 

Another useful quantity we shall consider is the profile, 
w;, of the interfacial adsorption, measuring the surplus of 
the non-boundary states in line, I, parallel to the ideal, 
straight interface. Then, 

W = ^ wi- 

i 


( 4 ) 
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Table 1. Numerical details used in the Wolff simulations of the 
8-state (r = 1/10) random-bond Potts model on the square 
lattice. In particular, the first column marks the linear size 
of the lattice, the second column the number of clusters TZd 
(created during averaging) and the third the number of the 
independent disorder realizations M. 


L 

Number of clusters TZc\ 

Disorder samples Af 

10 

2 X 10'^ 

20000 

20 

2 X 10'^ 

5000 

30 

2 X 10'^ 

4000 

40 

2 X 10'^ 

5000 

50 

30 X 10'^ 

4200 

60 

40 X 10'^ 

3200 

70 

50 X 10'^ 

1000 

80 

70 X 10'^ 

1000 

90 

200 X 10^ 

1000 

100 

300 X 10^ 

1000 


One may expect wi to fall off as one moves away from 
the interface. Its maximum indicates the location of the 
interface, which is, for a single realization of the Potts 
model with random couplings, not necessarily in the center 
of lattice. The spread, dw, of the adsorption profile may 
be measured by the distance, at which wi decreases to half 
of its maximal value, using a linear interpolation for the 
adsorption between successive lines, as will be discussed 
below. 

Using MC techniques, we recorded, in addition to the 
interfacial properties, standard thermodynamic quantities, 
for both types of boundary conditions. In particular, we 
measured the thermally averaged energy, Ei.,i and i?i: 2 , 
the specific heat given by the energy fluctuations, Ciu 
and Ci: 2 , and the order parameter given by the majority 
fraction of the Potts states mm, TOi:i and mi: 2 . 

In our simulations of the Potts models on square lat¬ 
tices with sites, we applied the Metropolis and the 
cluster-flip Wolff algorithm [^. Of course, cluster flips 
violating the boundary conditions are not allowed [53] . As 
usual, small lattices may be simulated using the Metropo¬ 
lis algorithm, while the Wolff algorithm is more efficient 
and is preferred for larger, say L > 30, system sizes. Over¬ 
all, we studied lattices with up to 100^ sites for the random 
Potts model and 200^ sites for its pure counterpart. 

Certainly, equilibration and averaging times depend on 
the lattice size. Moreover, for random models, we observed 
that the given bond realization may affect these times. In 
case of the Metropolis algorithm, eventually, simulations 
with 10^ Monte Carlo steps per site for L = 10 were per¬ 
formed, increasing the length of the runs, roughly, with 
L^. In case of the Wolff algorithm, the number of clus¬ 
ters T^ci used in our simulations varied from 2 x 10^ for 
the smaller systems sizes up to 3 x 10® for the larger sizes 
considered (see also Tabled]). 

The main source of error is, for random systems, the 
fact that the simulation data may vary quite drastically 
from bond configuration to bond configuration. The corre¬ 
sponding histograms or distributions have been recorded, 


especially, for the 8-state Potts model with r = 1/10. 
Bulk properties of this model have been studied quite ex¬ 
tensively before [19124125126127128] . Indeed, in the present 
study, we focused much attention on this case as well. The 
histograms at the critical point, for the various quantities 
discussed above, show nearly Gaussian shapes, but being 
weakly tailed, in accordance with previous observations 
and discussions for dilute Potts models on hierarchical 
lattices [5D]. The standard errors resulting from an en¬ 
semble average over bond realizations decrease with the 
number of configurations, Af, approximately, proportion¬ 
ally to Xj'/M. The proportionality factor seems to become 
somewhat smaller for larger lattices. To obtain reasonable 
accuracy, as will be elucidated below, we averaged over 
a large number of different bond configurations, as given 
explicitly in Table [T] For pure Potts models (r = 1) error 
bars follow from averaging over a few MC runs employing 
different random numbers, as usual. Typically, somewhat 
shorter runs are needed compared to those for the ran¬ 
dom models. Usually, the error bars are not depicted in 
the figures, because they are smaller than the sizes of the 
symbols. 

To determine critical properties from the MC data, 
we use finite-size scaling arguments. For example, for the 
interfacial adsorption, VF, one expects [3] 

IF « L“l7(tLi/^), (5) 

with a = \ — fijv and the reduced critical temperature 
t = \T — Tc\/Tc. 17 is the scaling function; /3 and v are 
the usual bulk critical exponents for the order parameter 
and the correlation length. A more refined ansatz invokes 
corrections to the asymptotic scaling behavior, as will be 
discussed in the following Section. 

3 Results 

3.1 Interfacial properties 

In the perfect 2D q-state Potts models the total interfacial 
adsorption W is known to vanish at zero and infinite tem¬ 
perature, with a maximum, for finite lattices, near the crit¬ 
ical temperature Tc- In the thermodynamic limit, W(T, L) 
diverges with characteristic critical exponents [318113] . In 
the present article, we shall take a closer look at the inter¬ 
facial adsorption by analyzing the corresponding profile, 
wi, as well. Moreover, the impact of randomness on the 
interfacial properties will be studied. 

Indeed, due to randomness in the couplings, the geom¬ 
etry of the interface between the I- and 2-domains may 
be rather complicated. Straightforward considerations on 
the ground-state energy show readily that, in contrast to 
the perfect case, the interface for a given random-bond 
configuration is at sufficiently small values of r, not nec¬ 
essarily straight or in the center of the lattice. Depending 
on that configuration, there may be even overhangs, i.e., 
interfaces may not satisfy the SOS criterion. As a conse¬ 
quence, the adsorption profile wi, for a fixed bond config¬ 
uration, may, at non-zero temperatures, display a maxi¬ 
mum away from the center of the lattice, as illustrated in 













’otts models on the square lattice 



Fig. 2. Critical adsorption profiles wi in the dilute 8-state 
Potts model, r = 1/10, for a single random-bond configuration 
(circles and squares) and for an average over N = 120 real¬ 
izations (triangles up). For comparison, the critical adsorption 

Ici in FViCi mi-rfi + "P/'iF'Fo io ^onin+tar^ /^Ffi o n rrlfio 



Fig. 3. Temperature dependence of the spread of the adsorp¬ 
tion profiles d„{T) in the perfect 3-state Potts model, as ob¬ 
served in simulations of lattices with 20^ (circles) and 60^ (di¬ 
amonds) sites. 


Figure [2] As shown in that figure as well, after averaging 
over the ensemble of bond realizations, wi is expected to 
be symmetric about the center of the lattice. 

The total interfacial adsorption W, at the critical point, 
for lattices of small and moderate sizes, tends to increase 
with the number of non-boundary states, q — 2, and it is 
smaller in the random case compared to the perfect model. 
This latter observation may be explained by the fact that 
in the random case the system is, both for predominantly 
strong and weak bonds, locally effectively below or above 
criticality, while the interfacial adsorption is expected to 
be maximal close to |3I8) . 

Figure [5] indicates that the spread in the adsorption 
profiles, dw, as defined in the previous Section, seems to 


increase fairly weakly due to randomness. It is found to 
depend rather mildly on q. Note that the temperature de¬ 
pendent spread dw(T) displays a maximum close to the 
critical point, as illustrated in Figure [31 similar to the 
behavior of W |3I8) . Its critical behavior, in the thermo¬ 
dynamic limit L —>■ oo, will be analyzed in the following 
subsection. 


3.2 Critical phenomena 

Following the finite-size scaling ansatz for continuous phase 
transitions, equation ®, the leading critical behavior of 
the interfacial adsorption, W, is given by [5] 

VF(Tc) oc ( 6 ) 

and 

IF(L —>■ oo,t) oc (7) 

with the critical exponents a and b being determined by 
the bulk critical exponents /3 and v [3] 

a = l —P/v] h = l3— V. (8) 

As stated before, the critical temperature, Tc, is known 
exactly for perfect and random Potts models on the square 
lattice, equation ©■ 

These predictions have been confirmed quite reason¬ 
ably in previous Monte Carlo simulations for perfect Potts 
models with q = 3 and 4. In the present study, we extend 
and refine the comparison by considering randomness as 
well and by including corrections to scaling. The (leading 
or effective) corrections may be cast in the form, as usual, 

IF(Tc) = IFoL“(l + coL-") (9) 

and 

W{L^CG,t)=Wit'^{l + city) ( 10 ) 

with t = |Tc — T|/Tc. In fact, the correction terms play an 
important role when fitting the MC data for the lattice 
sizes we simulated. 

Let us first discuss results for the perfect 3-state Potts 
model. Compared to previous simulations [3], the statis¬ 
tics has been improved significantly and larger lattice sizes 
have been studied. The size dependence of the critical ad¬ 
sorption, W{Tc), is depicted in Figure [H The simulation 
data for the complete lattice-size spectrum L = 10 — 200 
have been analyzed by fitting them to equation ([3]). We 
have also performed a fitting for the temperature depen¬ 
dence of the critical adsorption of the form (flUl) for the 
largest linear size studied, i.e., L = 200. The resulting esti¬ 
mates for the critical exponents, obtained from good qual¬ 
ity fittings with y^/dof G {0.6 — 1}, a = 0.867(2) and b = 
—0.728(6) agree nicely within errors with the predicted 
exact values a = 13/15 = 0.866 • • • (we remind the reader 
that P/iy = 2/15) and b = -13/18 = -0.722 •• • [IT]. 
Note that we fixed the leading correction-to-scaling expo¬ 
nents, X and y, to the predicted exact values, a; = 3/5 
and y = Aj3 |29I30I31] . Likewise, the estimates for the 
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Fig. 4. Finite-size scaling of the critical interfacial adsorption 
of the pure 3-states Potts model, showing the simulation data 
and the fitting curve. 



Fig. 5. Finite-size scaling of the critical interfacial adsorption 
data for the q = 8 and r = 1/10 diluted Potts model, showing 
the simulation data and the htting curve. 

exponents, x and y, are observed to agree with the pre¬ 
dicted values, when fixing the critical exponents a and b 
in the fittings. Accordingly, the findings on the perfect 3- 
state Potts model strongly support the correctness of the 
finite-size scaling description ([S]). 

The main focus of the current study is on the dilute 8 - 
state Potts model. Following previous considerations on its 
bulk critical properties, we set r = J 2 /J 1 = 1/10, where 
the randomness dominated behavior is expected to show 
up already for moderate lattice sizes 
In particular, we monitored, in our simulations, the size 
dependence of the critical interfacial adsorption W{Tc). 
As discussed above, the standard errors stem from av¬ 
eraging over many bond realizations. Numerical results 
are depicted in Figure [S] Fitting the MC data to equa¬ 
tion do]) with a correction-to-scaling exponent set to the 
value X = —1 which optimizes the fit, one obtains the 
value a = 0.856(7). 


Fig. 6. Size dependence of the critical spread of the adsorption 
profiles d„{T^) in the perfect 3- (circles) and 4- (squares) state 
as well as in the random, r = 1/10, 8- (triangles) state Potts 
models on square lattices with sites. Lines are simple guides 
to the eye. 


For the above random case {q = 8 and r = 1/10) and 
from previous MC simulations and transfer-matrix calcu¬ 
lations [19124125126127128) . fijv has been determined to be, 
approximately, 0.145 ±0.005. Actually, in the present sim¬ 
ulations, we also recorded the size dependence of the order 
parameter, for both types of fixed boundary conditions, 
TOi:i and mi: 2 , vanishing, at Tc, as ~ . The result¬ 

ing estimate for the exponent, based on simple power-law 
fits to the magnetization with systematically increasing 
the smallest lattice size, confirms the previous findings. 
Actually, a linear extrapolation of the fit exponent leads 
to the value fAjv = 0.145 ± 0.003 for toi:i as well as for 
toi: 2 . Accordingly, we may safely conclude for the random 
case as well a = 1 — (3/v, in accordance with the finite-size 
scaling ansatz, equation ©• 

We note in passing that our Monte Carlo data on the 
specific heat demonstrate that the absolute value of the 
corresponding critical exponent ajv is very small m- 
However, its more accurate determination is beyond the 
scope of the present study. 

Finally, critical properties of the adsorption profile, wi, 
and its spread, dw, will be reported. As indicated by Fig¬ 
ure 131 the spread dw at Tc grows when increasing the lat¬ 
tice size. We monitored the growth for the perfect 3- and 
4-state as well as for the dilute, with r — 1/10, 3-, 4-, and 
8 -state Potts models on the square lattice. The lattice size 
ranged from L = 6 up to L = 100. All data could be fitted 
well to the ansatz 


dw(dc, A) — do ± diL, (11) 


with the slope, di, depending rather strongly on the ran¬ 
domness, but only weakly on < 7 , as illustrated in Figure El 
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4 Summary 

We performed extensive Monte Carlo simulations to study 
critical interfacial properties in perfect and random ferro¬ 
magnetic q-state Potts models on the square lattice. In the 
dilute case, there are two distinct, strong and weak, bonds 
connecting neighboring sites. Both bonds are assumed to 
occur with the same probability, leading to self-duality. 
Interfaces have been introduced, by fixing the Potts vari¬ 
ables at opposite sites in two different states. The local 
Metropolis and the cluster-flip Wolff algorithms have been 
used. 

Randomness is found to affect, especially, the position 
of the interface, the excess or interfacial adsorption, and 
the form of the histograms resulting from the bond real¬ 
izations. 

Nevertheless, predictions of the isotropic finite-size scal¬ 
ing description for the interfacial adsorption at continuous 
phase transitions are observed to hold both for the per¬ 
fect and random cases, as has been exemplified for the 
pure 3-state and the dilute 8-state Potts models. In par¬ 
ticular, critical exponents of the interfacial adsorption are, 
indeed, determined by the bulk critical exponents for the 
order parameter and the correlation length. These bulk 
exponents depend on the number of Potts states, q. 

Moreover, we analyzed the adsorption profiles, describ¬ 
ing the spatial structure of the total interfacial adsorption. 
The spread of the profiles seems to diverge at the phase 
transition in the thermodynamic limit. For lattices with 

sites, we find a linear increase with the linear dimension 
of the lattice for perfect and random models, simulating 
the cases g = 3, 4, and 8. 
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